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\ INTRODUCTION 

For  years  the  advantages  of  being  able  to  produce  several  frequen- 
cies from  a single  master  source  have  been  recognized.  In  many 
applications,  the  easiest  method  of  producing  different  frequencies 
from  a master  clock  has  been  to  divide  the  clock  frequency  by  stages 
to  obtain  a frequency  which  is  a power-of-two  division  of  the  clock. 

This  method  has  the  intrinsic  advantages  of  extreme  simplicity,  stability 
commensurate  with  the  master  clock,  and  precise  repeatability.  Its 
disadvantage  is  that  only  a very  limited  number  of  such  frequencies  may 
be  produced. 

This  paper  describes  a method  by  which  a large  number  of  precise, 
stable  frequencies  may  be  derived  from  a single  frequency  source. 


SECTION  I 
DESCRIPTION 

Imagine  a register,  A,  of  length  N clocked  at  a frequency  folk 
(Fig.  1).  Imagine  a second  register,  F,  containing  a number,  M.  The 
A register  is  to  be  initially  set  to  zero.  Then  the  contents  of  the  F 
register  are  to  be  added  to  the  contents  of  the  A register  at  every  clock 
pulse,  and  the  result  is  to  be  stored  in  the  A register.  In  other  words, 
the  A register  is  to  be  run  as  an  accumulator.  The  rate  of  accumulation 
is  determined  by  the  number  M in  the  F register. 

Take  as  the  output,  the  contents  of  the  most  significant  Flip-Flop 
(FF),  Alt  of  the  A register.  Now  the  contents  of  the  FF  Ax  as  a function 
of  time  may  be  regarded  as  a square  wave  whose  frequency  is  deter- 
mined by  the  number  M carried  in  the  F register.  The  larger  the 
Output 
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The  F register  contains  a number,  M. 
FIG.  1 


number  M,  the  faster  the  rate  of  accumulation 


will  overflow,  and  the  faster  the  A 


It  will  be  clear  that  if  the  number  M is  a power  of  two,  then  the 
frequency  of  the  square  wave  output  will  be  the  clock  frequency  divided 
by  a power  of  two. 


The  A register  is  able  to  contain  the  number  2-1.  On  the  first 
clock  pulse  the  number  in  the  A register  changes  to  1.  On  the  2*  • 1 nl 
pulse  the  number  changes  to  2'  - 1.  On  the  2N  th  pulse,  the  A register 
overflows  and  the  contents  return  to  zero.  Thus  in  2*  pulses  the  Aj  FF 
has  gone  from  zero  to  one  and  back  to  zero  which  is  one  cycle.  The 
cycle  is  2*  pulses  long  and  the  frequency  is 


For  M set  equal  to  any  power  of  two  greater  than  zero,  say  21,  the 
A register  is  effectively  shortened  by  i Flip-Flops,  (Fig.  2),  as  all 
Flip-Flops  to  the  right  of  the  i**1  FF  never  contain  anything  but  zeros. 
The  output  frequency  is  now 


ARBITRARY  M 


Now  aisume  the  F register  to  contain  an  arbitrary  number,  M. 

Then  in  general  when  the  A register  overflows,  it  does  not  return  to 
zero,  but  to  some  positive  value  from  which  it  begins  a new  accumula- 
tion. 

Define 

l[x]  = the  integer  part  of  x 

G{x]  = the  fractional  part  of  x 

Example: 

x = 3.4  ^ I[xl  = 3;  C.[xl  = 0.4 

If  the  A register  begins  at  zero,  the  first  overflow  will  occur  in 
l^/M]  4 1 clock  pulses  and  will  leave  a remainder,  R,  in  the  A register. 
This  remainder,  R,  may  be  computed.  Note  that  the  number  M has  been 
added  to  the  A register  l[2N/M]  4 1 times.  In  other  words  an  attempt 
has  been  made  to  place  the  number  M(I[2N/Ml  4 1 ) in  the  A register. 

But  the  register  will  only  hold  the  number  2".  The  remainder  is  the 
difference: 

R = M(l{2"/Ml  4 1)  - 2*  (1-3) 

The  A register  will  overflow  in  l[2"/Ml  pulses  on  the  k 4 1st  cycle 
where 


kR  > M - R (1-4) 

and  will  continue  to  overflow  in  l[2"/M]  cycles  until  the  accumulated 
remainder  goes  to  zero,  i.  e.  , in  g cycles  where 


f 


g(M  - R)  = kR  (1-5) 

At  this  point,  the  A register  contains  the  number  zero  and  the  entire 
process  starts  over  again. 

Actually,  the  process  may  not  be  quite  so  simple.  The  remainder 
may  not  go  directly  to  zero,  and  the  A register  will  go  through  some 
sequence  of  overflows  in  l[2N/Ml  4 1 pulses  and  in  l[2N/Ml  pulses  until 
finally  the  remainder  does  go  to  zero  (and  eventually  it  must  go  to  zero), 
whereupon  the  entire  sequence  repeats. 

Each  time  the  overflow  occurs  in  l[2*/Ml  4 1 pulses, 

R = M(l[2"/Ml  4 1)  - 2N  (1-6) 

is  added  to  the  remainder.  Each  time  the  overflow  occurs  in  l[2N/M] 
pulses, 

M - R 

is  subtracted  from  the  remainder.  The  entire  sequence  repeats  when 
the  remainder  goes  to  zero,  as  that  is  the  initial  condition.  If  the  total 
number  of  times  the  overflow  occurred  in  l[2N/M]  4 1 pulses  is  given  by 
k,  and  the  total  number  of  times  the  overflow  occurred  in  l{2* /Ml  pulses 
is  given  by  g,  then  it  is  clear  that 

g(M  - R)  = kR  (1-7) 


as  before. 

Note  that  the  fraction  of  the  time  that  the  overflow  occurs  in 
l[2N/Ml  4 1 pulses  is  just 
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k _ M - R 
g + k M 


= g[zVm1 


(1-8) 


The  average  frequency  which  is  put  out  by  FF  Al  is  given  by 


f = 2n 


(1-9) 


The  average  frequency  is  arrived  at  by  the  output  switching  between  two 
frequencies,  namely 


“h  = l[2VMl 


(1-10) 


tl8* = ^TvmTTT 


d-ii) 


with  the  fraction  of  the  time,  t,  spent  at  the  latter  frequency  given  by 


t = G[2"/Ml  (Fig.  3) 


(1-12) 


The  separation  between  these  two  frequencies,  Af,  is  given  by 


I[2"/Mla  Vl[2N/M] 


(1-13) 


and  the  fractional  deviation  between  the  two  frequencies  is  given  by 


Af  2^ 

f.„  ‘ M(l[2"/Mla  + I[2n/M]) 


for  2N  » M. 


(1-14) 


(1-15) 
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l[2  /M]  4 1 


If  2 /Ml 


l[2  /M]2  + If 2 /M] 


FIG.  3 


LIMITATIONS 


The  output  waveform  may  be  thought  of  as  the  zero-crossings  of  a 
sine  wave  of  frequency 


f = MW 


which  is  then  sampled  at  the  clock  frequency.  Indeed,  this  method 
must  be  considered  as  a sampled  system.  It  should  be  clear  that  no 
frequency  greater  than  falk/2  may  be  obtained,  and  no  frequency  lower 
than  felk/2N  may  be  obtained  with  a register  of  length  N. 

The  precision  to  which  any  desired  frequency  may  be  produced  is 
determined  by  the  length  of  the  A register.  The  frequencies  which  may 
be  obtained  are  all  multiples  of  f el k /2N , so  that  any  frequency  may  be 
approached  to  within  fclk/2N+l  Hz.  A higher  precision  may  be  obtained 
by  lengthening  the  A register. 

In  theory,  any  frequency  below  folk/2  may  be  approached  to  within 
any  degree  of  accuracy,  simply  by  making  N sufficiently  large.  In 
practice,  however,  the  maximum  length  of  A,  and  therefore  is 

determined  by  the  circuitry  used  to  perform  the  addition  of  the  A and  F 
registers.  The  carry  produced  by  the  sum  of  the  least  significant 
digits  must  be  able  to  propagate  through  all  the  sums  to  the  most  signif- 
icant digits  during  the  time  interval  between  adjacent  clock  pulses, 

This  condition  sets  the  ratio  between  fclk  and  N,  and  effectively  sets  the 
upper  limit  on  the  range  of  frequencies  which  may  be  produced  by  this 
method. 

One  easy  and  simple  method  of  determining  the  approximate  fre- 
quency which  will  be  obtained  by  any  number  in  the  F register  is  demon- 
strated in  Fig.  4. 

Each  position  in  the  F register  corresponds  to  a frequency  which 
is  a power-of-two  division  of  the  clock.  A "1"  in  any  position  will  add 
the  frequency  corresponding  to  that  position  to  the  sum  of  the  frequencies 
corresponding  to  the  other  "l"s  in  the  F register. 
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The  output  frequency  is  given  by 


f = E -^5*  over  all  i's  for  which  a "1"  appears  in  the  F register 


Example: 


clock  = 1 mHz 
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1 

fa:  125  + 31  + 15.5  + 4 kHz  = 175.5  kHz 
f.„e,  = 175.781  kHz 


FIG.  4 


CONCLUSIONS 


The  great  usefulness  of  this  method  of  producing  square-wave 
frequencies  is  its  overwhelming  simplici'y.  The  frequencies  which  are 
produced  are  made  with  the  repeatability  and  stability  of  the  master 
clock,  and  may  be  identically  phased  simply  by  resetting  the  correspond- 
ing A registers  together. 


SECTION  II 


FREQUENCY  SPECTRA 

The  description  given  thus  far  indicates  that  the  output  waveform  is 
a square  wave  which  contains  frequency  components  that  are  not  present 
in  a "pure"  square  wave  of  the  desired  frequency.  The  equivalent 
circuit  for  the  system  is  shown  in  Fig.  5.  In  order  to  analyze  the  out- 
put spectrum  it  would  be  advantageous  at  this  time  to  derive  an  analytical 
expression  for  the  spectrum  of  the  output  waveform,  0(u.-). 

The  input  waveform  is  a "pure"  square  wave  of  the  desired  fre- 
quency, f(t).  This  function  of  time  has  a I ourier  transform  given  by 

J[f(t)l  s FU)  = j f(t)e"lundt  (2-1) 


FIG.  5 
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The  function  f(t)  is  sampled  at  times  T,  = 2 ^/td,  and  the  resulting  wave 
form  will  be  called  f*(t),  with  Fourier  transform  F*(u).  The  function 
f*(t)  is  fed  into  a holding  circuit  whose  output,  o(t),  with  Fourier  trans 
form  0(a:)i  is  the  function  which  the  frequency  generator  .actually 
produces. 


Given  an  input  function  f (t ) , we  have 


f (t)6 (t  - nT.)  = f(t) 


can  be  written  in  a Fourier  series  as 


It  follows  that 


and  its  Fourier  transform  is  indicated  by  the  formula 


so  that 


[f(t)e**]  = Fla  - a) 


?[f*(t)]  = F*M  = Y ^ F(<aj  - mw,)1 


(2-7) 


(2-8) 


Thus  it  is  shown  that  the  spectrum  of  f*(t)  is  just  the  entire  spec- 
trum of  l(t)  repeated  at  intervals  of  u:,. 

Given  the  following  sequence  to  compute  F(j^),  the  Fourier  trans- 
form of  f(t): 

if 


f(t)  = ^ ejBU,#t 


(2-9) 


where 


and 


*0  = 2ff /T 


1 ( '/J 

Of,  = ^ J_T/jf(t)e',,,u)oldt 


(2-10) 


then 


F(u;)  = 2ff  orB5(-a  - n^ 


(2-11) 


iLindorff,  David  P.  , "Theory  of  Sampled-Data  Control  Systems," 
Wiley,  N.  Y.  , 1965,  p.  30-31.  The  derivation  given  above  follows 
Lindorff's  text,  but  uses  the  notation  of  Papoulis  cited  below. 

^Papoulis,  Athanasios,  "The  Fourier  Integral  and  Its  Applications, 
McGraw-Hill,  N.  Y.  , p.  43. 
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The  actual  input  waveform  is 


-T/4 


T/4  3T/4 


Assume 


where 


then 


£(t)  = 1 - T/4  < t < T/4 


= 0 T/4  < t < 3T/4 


f(t)  = a.ejB~* 

« t — OB 


= 2rr/T 


*■  * 7 L 


i r '^4 

«■  * T /-/• 


-1 


jTnur0 


r */4 

e" 

L-t/4 


-4 


4jTnu.o 


1 sin  (ng,-aT/4) 

2 nu;0T/4 


(2-12) 

(2-13) 

(2-14) 

(2-15) 

(2-16) 

(2-17) 

(2-18) 
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Substituting  T = 2Wx0 


sin  (n?r/2] 


(2-19) 


and  its  Fourier  transform  is,  from  Eq,  2-11 


an 

F(w)  = 2ff  7 — sin  (n77/2)6(j.  - nag)- 

nrr 


(2-21) 


Thus  from  Eq.  2-8 

od  an 

F*(u)  = ~ ^ ^ ” sin  (nff/2)5U  - mA|  - nu.0)  (2-22) 

ms  -as  «r-  • 

Examining  the  clamp  circuit,  it  is  obvious  that  it  has  the  property 
that  for  an  impulse  input,  the  output  is  a square  pulse  of  a height  equal 
to  the  impulse  strength  and  of  width  T„#ll  = T,.  In  other  words,  the 
output  is  the  convolution  of  6(t)  and  a square  window  of 


= T#  = 2^/0,,. 


(2-23) 


The  Fourier  transform  of  such  a window  is  given  by 


-T./2 


and  combining  this  with  Eq.  2-13  yields 


an 

f(t)  = ^ — sin  {nn/2)ei'",*x 

Z-  nrr 


(2-20) 
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h(t)  = 1 - T,/2  < t < T,/2  (2-24) 

0 everywhere  else 


H(a) 

= f ' e',a,tdt 

i rT*/j 

= T~  e'4uH  / 

Jw  _T,/J 

(2-25) 

(2-26) 

_ zl  (e-ju.'T,/a  _ eJu}t,/8) 

(2-27) 

2j 

(2-28) 

= - sin  UT./2) 

uC 

(2-29) 

H(u>) 

sin  xtr/a. 

* a. 

(2-30) 

Now 

(2-31) 

(2-32) 

This  result,  while  interesting  as  a complete  expression  of  the  out- 
put spectrum,  is  somewhat  difficult  to  use  in  practice.  Since  the  input 
waveform  is  a square  wave,  it  contains  spectral  lines  all  the  way  out 
to  infinity.  It  is  clear  that  some  of  these  lines  will  be  folded  about  the 
line  0.(72  back  down  into  the  region  below  u',/2.  It  would  be  useful  to 
have  an  expression  which  bounds  the  fraction  of  the  average  power  which 


0(a)  = F*(a)H(a) 


0(a)  = 2tr  ~ sin  (n7T/2)<S(a  - ma,  - na0) 

u>  TT  / jjg  Lm*  Zl“ 
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is  folded  down,  and  also  one  which  bounds  the  most  prominent  spectral 
line  which  is  folded  back.  Such  an  expression  will  now  be  derived. 

For  any  function  of  time  which  is  expandable  in  a Fourier  series,  the 
average  power  in  the  waveform  is  given  by: 


f3(t)  * average  power  = N |o,  |2 

^ 4 


CD 

fa(t)  = Q!02  + 2 ^ | a,  |2  by  symmetry 


O',  = — sin  (n/r/2) 
n IT 


«■(*>■'«>*♦£  £ 3 

n odd 


2 , £ / , 1 1 1 V 

0,0  ff2  ( 1 + 32  + 52  + 72  + • • •) 


= Oo 


a + i.  ( l!)4 

1 tr2  ' 8 ' 


= Ofo3  + 1/4 


(2-33) 

(2-34) 

(2-35) 

(2-36) 

(2-37) 

(2-38) 

(2-39) 


The  average  power  in  that  portion  of  the  waveform  which  is  not 
folded  is  given  by  the  finite  summation: 


f3(t). 


1 


n odd 


noIflfluCo  < a.,/2 


(2-40) 


SPapoulis,  op  cit,  p.  241 

^Dwight,  Herbert  Bristol,  "Tables  of  Integrals  and  Other  Mathe- 
matical Data,  " Macmillan,  N.  Y.  , 1947,  p.  11,  Eq.  48.4. 


So  that  now  the  fraction  of  the  average  power  which  goes  into  the 
unwanted  spectrum  which  is  folded  back  is  given  by 


f3(t)u  = f*(t)  - fs(t). 


(2-41) 


"o 

= 1 /4  - — r for  n odd 

it*  na 


(2-42) 


But  this  last  summation  is  finite,  and  may  be  done  for  any  given  case 
in  a reasonable  length  of  time. 

The  fraction  of  the  power  which  is  folded  back  is  bounded  by  the 
ratio  of  the  unwanted  power  to  the  average  power: 


"ft 

— i 

qij). . * V n 

Pit)  ‘ 1/2 


JV.noj.0  < u.,/2 


(2-43) 


This  expression  is  an  upper  bound  of  the  power  which  goes  into  the 
folded  spectrum. 

The  most  prominant  line  in  the  folded  spectrum  is  just  the  lowest 
frequency  line  in  the  unwanted  spectrum,  namely  that  given  by  n = n«,  + 2. 
The  ratio  of  the  magnitude  of  this  line  to  the  magnitude  of  the  fundamen 
tal  of  the  square  wave  is  just 


FIG.  6 


which  produces  a frequency  of  69.336  kHz.  The  separation  between 
frequencies  is  0.9765  kHz. 

If  this  frequency  separation  were  not  sufficiently  large,  an  alternate 
approach  would  be  to  make  the  A register  a little  longer  and  space  the 
frequencies  a little  further  apart.  Set  N = 14.  Then  the  smallest  fr 
quency  steps  would  be  244.14  Hz.  The  smallest  M would  be  given  by 


30  kHz 
0.24414  kHz 


122.8 


(3-4) 


Therefore  set  M = 123,  which  corresponds  to  a frequency  of  30.029  kHz. 
The  next  frequency  would  be  31.250,  an  increase  of  1.221  kHz,  which  is 
produced  by  M = 128  (an  increase  of  5).  The  highest  frequency  would  be 
given  by 


M = 123  + 39(5)  = 318  f = 77.636  kHz  (3-5) 
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Using  Eq.  2-43,  that  percentage 

of  the  total  power  that  is  folded 

back  may  be  calculated.  This  gives 

n^n^gCo  >*.« 

u.o  = 77.636  kHz 

x,!2  - 2,000  kHz 

2,000 

n°  < 77  636  = 

: 25.761 

Therefore,  setting  n0  = 25 

no 

zL  na 

l 

is  given  by  the  following  table: 

n 

l/ns 

1 

1.00000 

3 

0. lllll 

5 

0.04000 

7 

0.02041 

9 

0.01234 

11 

0.00826 

13 

0.00592 

15 

0.00444 

17 

0.00346 

19 

0.00277 

21 

0.00227 

23 

0.00189 

25 

0.00160 

Vl/n3  = 

1.21447 

Thus,  Eq.  2-43  becomes 


1/4  - 2/ra  (1.21447) 
1/2 


0.5 


4. 85788 
9.86959 


= 0.5  - 0.49221 


= 0.00779 

= - 21  db 

And  the  most  prominant  line  is  down  by  1/27  compared  to  the  fundamen- 
tal, or  0.0370  = -28.6  db. 

Since  the  worst  case  is  represented  by  the  highest  frequency  made, 
the  calculation  above  sets  a maximum  upper  bound  on  the  TOTAL  folded 
power  and  the  most  prominant  spectral  line  for  all  the  frequencies  in 
this  example. 

The  author  used  Texas  Instruments  SN7483N  Quad-Adders  to  per- 
form the  addition.  These  components  have  a propagation  delay  in  the 
carry  term  of  12  nsec  maximum  per  addition.  The  time  between  clock 
pulses  (from  trailing  edge  to  leading  edge)  was  set  at  185  nsec.  Thus 
the  maximum  number  of  additions  which  ma;  be  safely  performed  is  15. 
This  determines  the  maximum  length  of  the  A register. 

A number  of  frequency  generators  utilizing  the  above  principles 
have  been  designed.  Each  frequency  generator  as  actually  constructed 
consists  of  12  dual-in-line  IO  chips  and  occupies  a volume  of 

V = 3 5/8  *2  1/8  * 1/4  inches3  = 1.9  inches3. 
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